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Abstract 

We reanalyze Oklo 149 Sm data using realistic models of the natural nuclear reactors. Disagree- 

i> . 

ments among recent Oklo determinations of the time evolution of a, the electromagnetic fine 
q_) , structure constant, are shown to be due to different reactor models, which led to different neu- 

tron spectra used in the calculations. We use known Oklo reactor epithermal spectral indices as 
criteria for selecting realistic reactor models. Two Oklo reactors, RZ2 and RZ10, were modeled 



with MCNP. The resulting neutron spectra were used to calculate the change in the 149 Sm effec- 



tive neutron capture cross section as a function of a possible shift in the energy of the 97.3-meV 
resonance. We independently deduce ancient 149 Sm effective cross sections, and use these values 
to set limits on the time-variation of a. Our study resolves a contradictory situation with previous 
Oklo o-results. Our suggested 2a bound on a possible time variation of a over two billion years is 
stringent: —0.24 x 10~ 7 < — < 0.11 x 10~ 7 , but model dependent in that it assumes only a has 
varied over time. 
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I. INTRODUCTION 

Two papers on the determination of the time evolution of a the electromagnetic 

fine structure constant from Oklo reactor data recently appeared, adding contradictory 
results to earlier investigations . The fractional change of a over a two billion year 

period has been found from Oklo data to be: Aa/a = +0.45ln'o7 x 1CT 7 [1], —0.56 x 1CT 7 < 
Aa/a < 0.66 x 10" 7 [2], -0.9 x 10" 7 < Aa/a < 1.2 x 10" 7 [sj] , and either -0.18 x 10~ 7 < 
Aa/a < 0.02 x 10" 7 or Aa/a = +0.88 ± 0.07 x 10" 7 in [6J. By comparison, astrophysics 
determinations from data on the shifts of the absorption lines in the spectra of quasar light 
have yielded Aa/a = — 72j^| x 10 -7 ?| and Aa/a = — 6lg x 10~ 7 [8] over an approximately 
ten billon year period. The sign of Aa is defined by the relationship Aa = a past — a pr esent, 
so that a negative sign, for example, means that 2-10 billion years ago the value of a was 
smaller than at present. For more results and references on the time variation of fundamental 
constants see Refs. 

As the results indicate, the situation is not entirely satisfactory: some analyses give only 
upper limits, while those showing a definite effect disagree even in sign. While theoretical 
models have been proposed which can accommodate time dependent rates of change of a, 
clarifying the disagreements among the Oklo analyses is important, particular since there 
are also questions about just how model-dependent these very precise limits actually are 
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131 ]. In this paper we will concentrate on the nuclear physics aspects of the Oklo 



reactors, focusing in particular on realistic models of the neutronics. 

The Oklo phenomenon has been known since 1972. The history of the discovery, the 
geological background, the relevant petrography, mineralogy, isotopic chemistry and the 
Oklo reactors physics are definitively described by Naudet 14| ■ Most of details of the Oklo 
phenomenon to which we will refer are from this largely unknown text. Findings from more 



recent Oklo studies are reported in Refs. 



15| and [16|. 



Sixteen natural uranium reactors have been identified in Gabon, West Equatorial Africa, 
in three different ore deposits: at Oklo, at Okelobondo 1.6 km away, and 20 km south of 
Oklo at the Bangombe. Collectively, these are called the Oklo Fossil Reactors. Well studied 
reactors include Zone two (RZ2) with more than sixty bore-holes, and more recently Zone 
ten (RZ10) with thirteen bore-holes. In RZ2, 1800 kg of 235 U underwent fission over 850 
kyr of operation and in RZ10 about 650 kg of 235 U fissioned (more rapidly) over 160 kyr 
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of operation. All reactor zones were found deficient in 235 U, and in most of them fission 
products were well retained. Isotopic concentrations were measured by mass spectrometry, 
and provided information on the neutron fluency, the neutron spectral index, and the 
235 tj restitution factor (burned 235 U is partly regenerated after a-decay of 239 Pu formed in 
neutron capture on 238 U). 



Due to the low 0.72%, abundance of 235 U and the high np capture cross section, present- 
day natural uranium cannot sustain a nuclear chain reaction with light water as a moderator. 
However, 2000 million years ago , when fission chain reactions started at Oklo, 235 U had a 
relative abundance of 3.7%, comparable to the 3—5% enrichment used in most commercial 
power reactors. In those times therefore a chain fission reaction was possible in principle and 
actually took place. Reactors in the northern part of the deposit, including RZ2 and RZ10, 
operated at a depth of several thousand meters, under then-marine sediments which came 
close to, but still below, the surface after the tectonic uprising about 250 millions years 
ago. At this depth, the conditions of pressure and temperature are close to those of the 
Pressurized Water Reactors (PWR) of today (temperature around 300 C, pressure about 20 
MPa). Of course, the Oklo reactor powers of 10—50 kW are greatly below the 1000-MW 
scale of the present reactors, and furthermore probably did not operate continuously. The 
authors of Ref. [17| deduced that RZ13 operated for a 0.5 hour until the accumulated heat 
boiled away the water, shutting down the cycle for up to 2.5 hours until the rocks cooled 
sufficiently to allow water saturation to initiate a new cycle. 

n 

Shlyakhter |3J was the first person to point out that a change in a could shift the position 
of the 97. 3- me V neutron resonance in 149 Sm and that as a result the present-day capture 
cross section could be different from the ancient value. Assuming a reactor temperature 
of 300K, and taking the fission isotope abundances known at that time, he found no 
evidence for a shift in the resonance position and accordingly obtained an upper bound 

or the fractional change in alpha of O.lxlO -7 (a revised number from comments in Ref. 

2I). Using updated abundance and temperature data, Damour and Dyson j^j], and later 
Fujii et al. [6| carried out more detailed studies for RZ2 and RZ10. They calculated the 



14( and 1.8 BY is used in Ref. 



1 The age of the Oklo natural reactors is debated in the literature: Fujii et al [6| and Dyson and Damour 
use 2 BY. Naudet cites 1.95 BY 
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present-day effective cross section by averaging the resonance cross section over a presumed 
fully thermalized Maxwellian neutron spectrum. In such an approach there is no need for 
a particular model for the Oklo reactors since the spectrum is determined solely by the 
temperature. Their results for the variation in a were basically in agreement, indicating 
no change. By contrast, in the recent papers [l| and [2j, where contradictory results have 
been obtained, the averaging is performed over neutron spectra with a 1/E epithermal 
tail in an addition to the Maxwellian contribution. Spectra with different contributions 
"rom the epithermal neutron tail were obtained with an infinite reactor model in Ref. 

i n 

JJ and from Monte Carlo modeling of a finite reactor in Ref. [2|. Not surprisingly, the 
use of different neutron spectra can lead to different results. But since these models are 
not unique, the question arises as to how to choose between them and between other models. 



In the present work we suggest using measured Oklo reactors epithermal spectral indices 
as criteria for selecting realistic reactor models. We perform MCNP calculations to find full 
scale models of RZ2 and RZ10 satisfying these criteria, and we use the resulting neutron 
flux spectra to calculate the dependence of the effective 149 Sm capture cross section on the 
resonance shift. We deduce independently the ancient 149 Sm effective neutron capture cross 
section using an updated formalism. From our limits on the 97.3-meV resonance shift, and 
assuming that only electroweak physics is varying, we we can set stringent limits on the 
time-variation of a from the Oklo data. 

The paper is organized as follows: In section [Til we review the definition of the effective 
cross section and the definitions of the various spectral indices used to define the contribu- 
tions of the epithermal neutrons to the neutron flux. In section 1111} we present our models 
of the Oklo reactor zones, and our MCNP calculations of the neutron spectra. In section 
ITVl we calculate the 149 Sm capture cross section as a function of the resonance energy shift, 
and in|V]we review our new calculation of the ancient 149 Sm cross section. In section PVTl we 
present our results and conclusions. 



II. EFFECTIVE CROSS SECTIONS AND SPECTRAL INDICES 

To analyze the Oklo reactor data without explicitly specifying the neutron den- 
sity n(v) at velocity v, it is customary to use not the average cross section a = 
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J °° n(v) a(v) vdv/ / °° n{y) vdv, but instead an effective cross section defined as: 

POO 

a = / n(v) a(v) v dv/nvo. (1) 



o 



Here n = n(v)dv is the total density, and ^o=2200 m/sec is the velocity of a neutron at 
thermal energy 0.0253 eV. To keep R, the reaction rate, unchanged it is necessary also to 
introduce an effective neutron flux density $ = / °° n(v) Vo dv different from the 'true' flux 
$ = J °° n(v) v dv, leading to R = a l> = a $. 

When the cross section for a particular reaction channel exactly follows a 1/v law 
(a = <7oVo/v), we have a = ao and the reaction rate <r$ = (Jquvq does not depend on 
the temperature T. When the cross section a deviates from the 1/v law at low energies (as 
it does for 149 Sm), and when the neutron spectrum is not pure Maxwellian (as is the case 
in any realistic reactor), the effective cross section can be written 

a = g{T)a + r Q I, (2) 

where g(T), a function of the temperature T, is a measure of the departure of o from the 
1/v law, J is a quantity related to the resonance integral of the cross section and To is the 
Oklo reactor spectral index, a measure of the contribution of epithermal neutrons to the 
cross section. The more well known epithermal Westcott index r [3| is related to tq by 

ro = r ^¥ (3) 

The Westcott index is a temperature dependent quantity while ro (as shown below) and / 
are independent of temperature. 

Since we will be concerned with extracting rp from neutron spectra calculated by MCNP 
for specific reactor models, we follow Ref. [18] and introduce the total neutron density n(v) 
and its epithermal fraction f ep i. 

n(v) = n [(1 - f epi ) n th (v) + f epi n epi (v)] , (4) 

where n t h(v) = (4/7r)(v 3 /v^) exp(— v 2 /v\) is the thermal Maxwellian distribution, and 
n ep i{v) = v c /v 2 for v > v c , otherwise zero, is an epithermal distribution which holds for 
systems with zero resonance absorption. Each of the distributions n t h{v) and n epi (v) is nor- 
malized separately to unity. The velocity v c is an as-yet unspecified cutoff velocity for the 
epithermal distribution and Vt is the most probable neutron velocity for Maxwellian with 



the temperature T as given by = v (T/T )2. The temperature T is the temperature of 
the Maxwellian density distribution having the most probable velocity vq, defined from the 
relation mv%/2 = kT = 0.0253 eV. 

The authors of Ref. [l^] suggest cutting off the epithermal distribution at energy E c = 
fikT with the parameter value \i ~ 5. The corresponding temperature dependent threshold 
velocity is then v c = Vt-JJL- At such a threshold the Maxwellian component is already 
several times larger than the epithermal one, so such an approximation is satisfactory. The 
Westcott spectral index r is then defined as 



fepi(T^) — f e pi- (5) 

To assess the temperature dependence of f epi (T), we transform from neutron densities to 
neutron fluxes. This yields for the total thermal component <p t h = n(l — f ep i)v = n(l — 



fepi) v T\J and for an ideal (without resonance absorption) epithermal flux per unit of 
neutron lethargy (lethargy is u = ln(10/-EMev)) 0epi(Aw = 1) = nf epi v c /2. Introducing the 
ratio 

^ = (frepijAu = 1) ^ 

4>th 

we find: 

^ y/^f^ fepi /j^ 

4 1 fepi 

For small values of f epi and /i ~ 5, we see that 5 ~ f epi . The quantity 5 is calculated in 
reactor physics [18] and for the case of moderation by hydrogen is 



6= m)i^r- (8) 

Here S a (u ) is the summed macroscopic absorption cross section at neutron velocity t>o, 
and S(-ff) is the macroscopic scattering cross section of hydrogen in the epithermal region. 
From this equation and identifying 5 with f epi we see that the epithermal fraction of the 



neutron density in a reactor behaves as f ep i(T) oc yT /T, thereby confirming that the Oklo 
spectral index defined by Eq. [3] is independent of temperature. 

In realistic systems the epithermal flux deviates from the 1/E dependence due to ab- 
sorption of neutrons in uranium resonances above an energy of several eV. In this case, two 
different definitions of ro, both approximations, have been used. Ref. [ljj] relies on the 
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shape of the neutron flux and defines 



ro 



'ept 



U 



(9) 



J th 



where <f) e pi{u*) is the is the flux per unit of lethargy at some effective energy in the resonance 
region and (f>* h is the total flux integrated up to an energy where the flux begins to increase 
above the 1/E level. In Ref. 18] the spectral index r at room temperature is computed as 
the ratio of the effective macroscopic absorption cross section S a to the moderating power 
££ s [18], and therefore 

ro = ^- (10) 

We compare these three definitions in the next section. In Ref. [1], the tempera- 
ture dependent Weinberg-Wigner 19( thermalization parameter A(T) = 2AY, a (kT)/'E s 
is used. For the case of hydrogen moderation only, it is related to ro by ro ~ Jn/AA(T )/2. 



The Oklo spectral indices ro are known quantities; their values have been deduced for 
several Oklo reactor zones [jjj], [3], [l(|, [20] from analysis of the fission products 143 Nd, 
147 Sm, and the 235 U concentrations. We will use these spectral indices to discriminate 
between possible models of the ancient reactor zones. In particular, the following exper- 
imental values have been deduced: ro = 0.20 — 0.25 for RZ2 and ro = 0.15 ± 0.02 for 
RZ10. The RZ2 result is a range of the borehole SC36 values corrected in Ref. 14| after 



re-evaluation of the cross section <3"i43 Nd to the expression a = 335 — 100r o . We report the 
RZ10 result as an average from four samples of Ref. [if] ]. 



III. OKLO REACTOR MODELS AND NEUTRON SPECTRA 

The reactor criticality is determined by the size, geometry and composition of the 
active zone. The latter influences the energy dependence of the neutron flux. The Oklo 
reactor zones include uraninite UO2, gangue (oxides of different metals with water of 
crystallization), water, and poisons that are present initially or build up during operation. 
Among these parameters, the most uncertain is the amount of water present at the 
time of reactor operation. In our modeling, we vary this parameter to match the experi- 
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mental spectral indices for RZ2 and RZ10 while keeping the reactor under critical conditions. 



The Oklo reactor cores have a characteristic horizontal size of order 10 m and occur 
in sandstone as lens shaped bodies of thickness varying between 20 and about 90 cm. 
The uranium content ranged from 20 to 80 wt.%. Each core is surrounded by a clay 
mantle. In the past water filled spaces left by cracks and fissures. The effective poros- 
ity required to achieve criticality is large (~20%), and is explained by a de-silicication 
process Ijj] consisting of partial leaching of the silica by hot thermal underground water. 



TABLE I: Composition (in g/cm 3 ) and neutronic parameters of the Oklo reactors. The values 
shown values are for 2 BY ago, for example, 30 wt. % of UO2 in RZ10 dry ore then corresponds 
to ~ 22 wt. % in present day ore. 



Zone 


uo 2 


H 2 


Si0 2 


FeO 


A1 2 3 


MgO 


MnO 


K 2 (Ca)0 


Total 


k eff , at 300K 


r f 


P 


RZ2 
RZ10 


2.500 
0.850 


0.636 
0.355 


0.359 
0.760 


0.149 
0.320 


0.238 
0.510 


0.077 
0.160 


0.009 
0.020 


0.020 
0.040 


3.99 
2.96 


1.033 * 
1.036 


0.22 
0.15 


0.800 
0.845 



' The ro values are calculated from the neutron densities (see text). 

* The value for the composition with a poison of 10 ppm Boron-10 equivalent. 



We obtained the Oklo reactor neutron fluxes at several temperatures with the code 



MCNP4C 



25] using the free gas option for the neutron scattering. Our model of an Oklo 



reactor is a flat cylinder of 70 cm height, 6 m diameter, surrounded by a 1 m thick reflector 
consisting of water saturated sandstone. The compositions (in g/cm 3 ) of the two reactors 
are shown in Table [fl Reactor zone RZ2 parameters are from Ref. 14J and reactor zone 
RZ10 components are from Refs. 15j and 16(. The total density of the active core material 
at ancient times was about 4 g/cm 3 for RZ2 and 3 g/cm 3 for RZ10. The most striking 
difference between RZ2 and RZ10 is the small amount of uraninite U0 2 in RZ10. In fact, 
with such a small amount of U0 2 it was not possible to make RZ10 critical with a poison 
more than 1 ppm of Boron-10 equivalent. The amount of water shown in Tabled] is a total 
amount, also including water of crystallization. 

To check the effects of chemical bonding, we repeated MCNP calculations using the 
option that considers scattering from hydrogen bound in water molecules. The spectra were 
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essentially unchanged except at the very lowest energies, and calculations of &uq (see section 
HVj) showed 1% changes or less at all temperatures. All calculations were therefore made 
with the free gas option. 



RZ10 MCNP fluxes and densities 




Neutron energy (eV) 



FIG. 1: Reactor zone RZ10 MCNP neutron fluxes and neutron densities for temperatures 20, 100, 
200, 300, 420, 500 C. The neutron fluxes are plotted as cpE and are the family of curves starting 
from a temperature of 20 C at the lower left. The neutron densities (normalized to one neutron 
per unit volume) are the families of curves starting upper left. 

The MCNP neutron fluxes and neutron densities are shown in figs [1] and [2] for temper- 
atures 20, 100, 200, 300, 420, 500 C. The neutron fluxes are the family of curves starting 
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at the lower left of each figure. The leftmost curve corresponds to a temperature of 20 C 
and the rightmost curve to 500 C. The uranium absorption resonances are prominent in the 
epithermal region. Also of note is that the flux is not flat in the epithermal region, indicating 
the spectrum is not precisely 1/E. The neutron densities (normalized to one neutron per 
unit volume) are the families of curves starting upper left in each figure. The topmost curve 
corresponds to 20 C and the lowest curve corresponds to 500 C. 

Oklo spectral indices can be calculated from these density plots, using Eqs. [3] and [5], with 
a cut off parameter \i — 5. The values we get for tq are shown in Table [H They agree with 
the experimental values cited earlier, confirming that we have realistic models of the reactor 
zones. The resonance escape probability p shown in the last column of Table fl] is discussed 
in Section [Vj 



TABLE II: Moderation characteristics of the Oklo reactor zone RZ2. For this reactor the hydrogen 



to uranium atomic ratio is ^f 1 - = 7.6, and the spectral index from Eq. [10] is found to be tq = 0.24. 



N v 



Atom 


Ni, 10 21 cm" 3 


Sa^cm" 1 




6 




235jJ 


0.207 


0.140 


0.003 


0.008 


0.0000 


238jJ 


5.374 


0.014 


0.048 


0.0008 


0.0003 


io b 


0.011 


0.042 


0.000 


0.210 


0.0000 


H 


42.60 


0.014 


0.873 


1.000 


0.7591 


O 


47.50 


0.0000 


0.204 


0.120 


0.0213 


Si 


3.57 


0.0006 


0.008 


0.011 


0.0001 


Fe 


2.92 


0.0030 


0.003 


0.036 


0.0001 


Al 


7.30 


0.0007 


0.005 


0.070 


0.0004 


Mg 


1.35 


0.0001 


0.005 


0.080 


0.0004 


Mn 


0.08 


0.0010 


0.0003 


0.036 


0.0000 


K 


0.26 


0.0005 


0.0005 


0.040 


0.0000 


TOTALS 


N = 111.2 


S a = 0.216 


S s = 1.150 




i = 0.782 



A specific feature of our present model as compared with Ref [1] is the presence of 
gangue in the reactor core. This means the reactor moderator is a composite substance with 
atoms besides hydrogen participating in the slowing down and absorption of neutrons. The 
contributions of these other atoms are presented in Table [III for RZ2, and in Table II III for 
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FIG. 2: Reactor zone RZ2 MCNP neutron fluxes and neutron densities for temperatures 20, 100, 
200, 300, 420, 500 C. The neutron fluxes are plotted as cpE and are the family of curves starting 
from a temperature of 20 C at the lower left. The neutron densities (normalized to one neutron 
per unit volume) are the families of curves starting upper left. 

RZ10. While only oxygen adds an appreciable amount to the moderating power ££ s , (here 
£ is the effective logarithmic energy loss defined as £ = J2i [£iE S i/£ s ]) several other elements 
contribute to the absorption parameter S a . Using these parameters, alternative values for 
ro can be calculated from Eq. [TOj as shown in the tables. These values agree well with the 
values determined from the neutron density spectra as shown in Table [H We confirmed also 
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that using the neutron flux spectra with Eq. [9j good agreement for ro values is found if the 
value for the flux per unit of lethargy is taken at about 100-eV neutron energy; this takes 
into account the deviation seen in the MCNP simulations of the flux from the 1/E law. 

TABLE III: Moderation characteristics of the Oklo reactor zone RZ10. For this reactor the hy- 
drogen to uranium atomic ratio is ^ = 13.0, and the spectral index from Eq. [TO] is found to be 
r Q = 0.14. 



Atom 


ivr. irv21 — 3 

iNj, iu cm 


^ -l 

^ai j Cm 


v -1 
2-*si, cm 


c. 

Ki 




235tt 
U 


U.UDo 


U.U4DU 


U.UU1U 


n none 
U.UUUo 


U.UUUU 


238tt 
yj 


1 77 


0048 


01 57 


0008 


0002 


H 


23.76 


0.0078 


0.4871 


1.000 


0.6667 





44.39 


0.0000 


0.1910 


0.120 


0.0314 


Si 


5.60 


0.0010 


0.0123 


0.011 


0.0002 


Fe 


2.94 


0.0076 


0.0076 


0.036 


0.0004 


Al 


7.29 


0.0017 


0.0017 


0.070 


0.0002 


Mg 


3.16 


0.0002 


0.0126 


0.080 


0.0015 


Mn 


0.18 


0.0024 


0.0004 


0.036 


0.0000 


Ca 


0.43 


0.0002 


0.0012 


0.040 


0.0001 


TOTALS 


N = 89.59 


S a = 0.072 


S s = 0.731 




£ = 0.701 



The calculated spectral index for the infinite medium reactor model of Ref. PJ is 
tq = 0.53, which is far away from the experimental value of 0.15. The difference is mainly 
due to the small amount of water (-^ ~ 3) which is much less than the values we found for 
RZ2, and RZ10. The possibility of a larger amount of water in the core was excluded on the 
grounds that the reactor would became divergent. While this is true for an infinite reactor 
including uranium, water and a poison, it is not true for a finite Oklo reactor that contains 
an additional constituent — a gangue. Our MCNP modeling confirmed that leakage of 
neutrons from a finite reactor of about 70 cm thick composed of only uranium and a small 
amount of water was greater than from an identical geometry reactor with gangue and 
much more water. As a consequence the reactor of the Ref. |lj] is undercritical if it is made 
finite. 
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The reactor models in Ref. Q| are close to ours in geometry but the active core 
compositions differ in uranium and water content. Spectral indices ro were not reported. 
The 'case 1' model of Ref. [2] has 1.5 g/cm 3 of UO2 and 0.355 g/cm 3 H 2 at a total density 
of 3.4 g/cm 3 with atomic ratio jj£ = 7.0. This ratio is below our value of 13.0, therefore 
we believe it does not accurately represent RZ10. The j^ 1 - = 7.0 ratio is close to the ex- 
pected value for RZ2 but the other compositions deviate from values established in Ref. [ijj . 



IV. THE CROSS SECTION AVERAGING PROCEDURE 

Following previous work, we evaluate the effective capture cross section a for 149 Sm 
numerically. The calculations were carried out with the code SPEAKEASY {22]. We include 
all resonances up to 51.6 eV along with the sub-threshold resonance at —0.285 eV. We neglect 



resonance interference terms, and use parameters from Mughabghab et al [23]. In practice, 
the contribution of the 97.3 meV resonance dominates; other resonances contribute only a 
few per cent to the sum. The bulk of the calculations were carried out neglecting Doppler 



broadening; separate calculations with the code SAMMY [24| indicated Doppler effects on 
(Tug were essentially negligible for all temperatures considered, as also found by earlier work. 

The MCNP calculations provide the neutron flux per energy bin, on a lethargy grid u 
running from 23 to 9.3 in steps of 0.1. This gives neutron energies E = 10 7 e~" in eV from 
about 1 meV up to 1 keV, with bin widths AE = 0.1E. A finer grid is required for the 
numerical integrations of cross sections for resonance neutrons since typical total resonance 
widths in 149 Sm are of order 100 meV. Accordingly, the MCNP flux is interpolated on to 
a lethargy grid with step size 0.001 (AE = 0.001.E). The flux is also re-normalized by 
requiring the total neutron density to sum to unity (n = 1). 
The effective cross section a is evaluated numerically using 

* = — (11) 

Here is the Breit Wigner cross section for resonance j at neutron energy E^, (pi is the 
neutron flux at energy Ei in a bin of width 0.001-Ej, and Vq = 2200 m/s is the velocity of a 
thermal neutron. 

The Breit Wigner cross section for a neutron of mass m, energy E^ at resonance j 
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Sm RZ2(dotted) and RZlO(solid) 




-0.2 -0.1 0.1 0.2 

Resonance shift (eV) 

FIG. 3: Calculations of the 149 Sm effective capture cross section 1T149 as a function of a possible 
resonance energy shift A r from -200 meV to +200 meV. The results shown are for RZ2 (dotted 
lines) and RZ10 (solid lines), and for temperatures T = 20, 100, 200, 300, 420, 500 C starting from 
the top. 

having resonance energy E r j, neutron width T ni j, and decay width T 7 j is given by: 

_ 90j 1T ^ 1 ^nij^-yj 

aij ~ 2mE i (Ei — E r j) 2 + r? /4' ^ 
14 



The energy dependent neutron width (assumed an s-wave) for resonance j is given by 



T nij = ^Ei(eV)T n0j} (13) 

where the product of the spin statistics factor g Q j and the reduced neutron width r ra0 j is 
tabulated in Ref. 



231 ] . The total width is given by 

Tij = T^j + T n ij (14) 

and the resonance energy is given by 

E r j = E r jo + A r (15) 

where E r j Q is the (present day) resonance energy, and A r is the energy shift associated 
with a possible change in the value of the fine structure constant a. 

Summing over all resonances, we evaluate 0149 for values of A r from -200 meV to +200 
meV. The results are shown in figure [3] for RZ2 (dotted lines) and RZ10 (solid lines) for 
temperatures T = 20, 100, 200, 300, 420, 500 C. Temperatures run from top to bottom; 
in each case T = 20 C is the top curve, and 500 C is the bottom curve. All resonances 
up to 51.6 eV are included in calculations, and are assumed to shift the same amount. 
Calculations with Maxwell-Boltzman spectra were also performed; in all cases they agreed 
closely with previous work. 



V. CALCULATION OF THE ANCIENT CAPTURE CROSS SECTION 

To rederive the ancient 149 Sm capture cross section, we consider the time evolution of 
the number densities A0i(t) of the six isotopes of interest - 235 U, 238 U, 239 Pu, 147 Sm, 148 Sm, 
and 149 Sm. We model the time dependence of the densities by the following set of coupled 
differential equations: 

dN " ~a 5 , tot N 5 + ^N 9 (16) 



(fidt 
dNs 
4>dt 

dN g „_ 0Ar A 9 
4>dt 



a° 8 N 8 - T (17) 
a° 8 N 8 - ^N 9 - a 9 , tot N 9 + T (18) 
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diVi47 



4>dt 
dN 148 



r = (i - P )(u 9 a ft9 N 9 + u 5 a ft5 N 5 ) (19) 

-0-147^147 + h,fY 5>U7 N 5 + ^9,fYg M7 Ng (20) 



-& 148^148 + & 147^147 (21) 



(j)dt 

-(TMgiVug + (7148^148 + a 5 jY 5;14 gN 5 + & 9 jY 9M9 N 9 . (22) 



4>dt 

In these equations, the subscripts 5, 8, and 9 refer to 235 U, 238 U, and 239 Pu, respectively, 
and tot and / refer to the total (absorption) and fission cross sections. The average neutron 
flux, 0, is assumed to be constant on a time scale long compared to the decay constants 
in the problem. The Y represent the fission yields to the Sm isotopes of interest, the 
subscripts 147, 148, and 149 refer to 147 Sm, 148 Sm, and 149 Sm, respectively, and u 5 ^ 9 are the 
average number of neutrons produced by 235 U and 239 Pu, respectively. We take the reactor 
start time as two billion years ago at to = and the end of operation for RZ10 as t\ = 160 
kyr and for RZ2 as t\ = 850 kyr. 

Our calculation differs from the work of Fujii et al. 6J in that we explicitly incorporate 
the resonance escape probability, p in modeling the restitution of 235 U. We also include a 
term for the 239 Pu alpha decay rate, l/Ag = 24.141/ ln(2) = 34.8 kyr which leads to a 
reduction in the restitution of 235 U, the issue being that 34.8 kyr is not small compared to 
the reactor operating lifetime of 160 kyr for RZ10. 

The resonance escape probability is the probability that a neutron produced by fission is 
not absorbed on the principal absorber in the system, 238 U, when moderating through the 
resonance region. The quantity (1 — p) therefore gives the fraction of fission neutrons that 
are absorbed by 238 U, most of which eventually converted to 239 Pu. The conversion factor 
C is defined as the ratio of atoms 239 Pu produced to the number of atoms 235 U "burned" in 
thermal neutron absorption. It is calculated (Eq. (7.8) in Ref. [3|) as 

C = RT 1 J*- + (l-p)v fi P*-, (23) 

&5,tot &5,tat 

where the first term is the contribution to conversion from thermal neutrons (here R is the 
fraction of atoms 235 U in the uranium fuel), and the second term is due to the absorption 
of neutrons by 238 U in the resonance region. Not all converted atoms end as the restituted 
235 U atoms, because some atoms 239 Pu will be burned in the thermal neutron flux <p. The 
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restitution factor C* is defined as 

C* = £ . (24) 

1 + 0^ 

The expected value of p is shown in Table [J and is calculated as follows. For U in a 
composite moderator, the definition of the resonance escape probability is 



1 r E„ dE" 



N U T res 



(25) 



where the effective resonance integral l r e e ^ for a homogeneous mixture of fuel and moderator 
was found in experimental and theoretical works [l9] (Eq. (10.29)) to be 

/1Q24V \ - 42 

i: "' =3 - 8 {^r) • (26) 

In these equations, S s is the total macroscopic scattering cross section of all elements in 
he active core, and £ is the corresponding effective logarithmic energy loss defined by (Ref . 
3,Eq. (10.20)) as £ = 

Note that these equations assume that the atoms are free for the energy loss, i.e., 
chemical bond effects are neglected. Because the resonance integral is over the range 
0.5 eV to 100 keV, the effects of chemical bonds are expected to be small. From these 
equations and using the RZ10 characteristics from Table IHIl we calculate for RZ10 the 
value (1 - p) = 0.155. Eq. El gives C(RZ10) = 0.11 + 2.05(1 - p) = 0.43. For the 
metasample RZ10 flux (see later) we then get C* = 0.77C = 0.33, in reasonable agreement 
with C* = 0.38, the average of the values quoted in Hidaka and Holliger (HH) jl^ . 



The coupled equations are solved with input parameters l_4j listed in (Table IIVI) . We 
confirmed the correctness of these effective cross sections by also calculating the integral a = 



J a(E)®(E, T)dE/ J $(E,T) /^jE/0.0253dE) with our MCNP fluxes. Because the resonance 
escape probability represents the capture probability for E > 0.5 eV, o"g is determined by 
integrating to 0.5 eV, the rest of the resonance integral being incorporated into the resonance 
escape probability. 

The coupled equations are solved as follows. First, fixing 1—p = 0.155 for RZ10, and with 
starting parameters iV 5 (0) = 0.0370, iVg(0) = 0.963, A^(0) = 0, we can solve the first three 
equations for the fluence ifiti which reproduces the N 5 (ti) / N 8 (ti) = (N 5 (ti) + N 9 (ti)) / N 8 (ti) 
ratio determined by the present-day HH data for each RZ10 sample. We use T Q = 2 x 10 9 
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TABLE IV: Effective cross sections for RZ10 and RZ2 used to determine CT149 (Cross sections are 
based on data from Naudet. Numbers in the first row are for the zone RZ10 and in the second row 



- for RZ2). 



0-5, kb 


<7 5 j, kb 


09, kb 


agj, kb 


a 8 °, kb 


(T147, kb 


<Ti4 8 , kb 


^5,147 


^9,147 


^5,149 


^9,149 


^5 


^9 


0.656 


0.549 


2.05 


1.35 


0.0027 


0.142 


0.0024 


0.0226 


0.0226 


0.011 


0.011 


2.43 


2.88 


0.665 


0.550 


2.16 


1.40 


0.0027 


0.184 


0.0024 


0.0226 


0.0226 


0.011 


0.011 


2.43 


2.88 



years ago to derive the ratios at time t\. We also calculate values for a metasample which 
represents an average of the individual sample properties. The idea of the metasample is 
that there could have been internal mixing between the samples or other processes. It is 
also useful to see whether the sample averages provide results in a reasonable range. 

With the fluence fixed, we next vary the starting amount of elemental Sm relative to 
U, -/Vs m (0)/iVu(0) to reproduce the ending isotopic fraction /W^i), and then vary the 
cross section cr 149 to reproduce the ending isotopic fraction /i4g(ti). We use the natu- 
ral abundancies to convert from an elemental starting ratio to the isotopic starting ratios 
(/ 144 (0) = 0.031 Ji 47 (0) = 0.151, /i4s(0) = 0.113, and /i 49 (0) = 0.139). In principle the 
starting fraction of Sm relative to U, while very small, is not a free parameter; it is deter- 
mined by data for /i 44 (ti) since 144 Sm is not produced in fission: 

fiu(ti)N Sm {h) = /i 44 (0)iVs m (0). (27) 

In practice, while <j 14 g can be tuned closely to match /i 49 (ii) we were not able to reproduce 
/i47(£i) without also varying iVs m (0)/jVu(0). Table IVl summarizes our results. The measured 
values from HH are shown above the empty line and our calculated values for N5(ti)/Ns(ti), 
/i47(ti) and /i4 9 (ti) are shown below the empty line, along with our calculated fluences, 
calculated starting fractions of elemental Sm relative to U, and calculated <5"i4g values for 
the four RZ10 samples and the metasample. The mean of the four values is 89.4 kb, with 
sample standard deviation 6.8 kb. 

The starting fractions, iVs m (0)/iVu(0) vary, but are consistent with HH data for bore 
holes 1640 and 1700 outside the reactor zones, and apart from 1480 are within a factor 
of two of values predicted by the HH data. The mean value for the ending amount of 
-^Sm(^i) = 74.6 ± 19.6 /igms/gm agrees well with the mean of the measured HH values 
67.8 ± 36.0 /igms/gm, also giving us confidence that the starting ratios are reasonable. 
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As a final value for RZ10 we adopt the cross section extracted for the metasample, and 
to take account of sample-to-sample variations we adopt the sample standard deviation as 
a la uncertainty: 

6-149 = 85.0 ± 6.8kb. (28) 

This is in good agreement with the mean of the four RZ10 samples analyzed by Fujii et al. 
jfj]: 91.2 ± 7.6 kb, again taking their sample standard deviation as a la error. 

Data for fifteen RZ2 samples were analyzed by Damour and Dyson [5], who set a 2a bound 
57 < 6"i49 < 93 kb. Results from RZ2 have tended to show more scatter because the samples 
come from mining near the surface with a greater potential for contamination compared to 
the deep underground samples for RZ10. Nevertheless we repeated our calculations for the 



five RZ2 samples from borehole SC36 - 1408 to 1418 - cited by Ruffenach 21| as having the 
most important contributions of elements formed in fission. We followed the same procedure 
as for RZ10, also including analysis of a metasample formed by averaging the input data for 
the five individual samples. With (1 — p) fixed, we varied the flux to match the U isotope 
ratios, and then varied 03,49 and the starting elemental ratio of Sm and U to reproduce the 
ending fx^{t\) and fimiti) fractions. Table IVT1 summarizes the results. The average of the 
five values is 75.6 kb with a sample standard deviation of 10.0 kb. The metasample value is 
71.5 kb, and we therefore adopt as our RZ2 result: 

0i49 = 71.5 ± lO.Okb, (29) 
a value consistent with the RZ10 result and with the Damour and Dyson analysis. 

VI. RESULTS FOR RESONANCE SHIFT AND CHANGE IN a 

To calculate bounds on the possible shift of the 149 Sm resonances over time we use the 
RZ10 result 0149 = 85.0 ± 6.8 kb, and conservatively adopt a ±20 range to establish limits: 
71.4 < 149 < 98.6 kb. We further assume that the reactor operated between 200 and 300 C 



as per the analysis by Meshick et al. [17J], which limits the temperature to this range, and 
as per knowledge that the geological formation and depth would not allow liquid water to 
be present above 300 C. 

An expanded version of the calculated RZ10 0149 cross sections at 200 and 300 C is shown 
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TABLE V: RZ10 calculated fluences <j)t\ and ancient cross sections <7i4g. Input data for Sm and U 
isotope fractions (above the empty line) are from Ref. [if]]. Calculated values are below the empty 
line. For these calculations t\ = 160 kyr. The metasample input data are an average of the data 
for the four individual samples. 



Property 


SF84-1469 


SF84-1480 


SF84-1485 


SF84-1492 


Meta 


N 5 (ti) 
N a (ti) 

7144^1 ) 


U.Uol i D 

n nrn n^9 


U.UZDOZ 


H H9Qfi7 

u.uzyo / 

009073 
U.UUiU / o 




u.uzyoz 

U.UUl 1 ou 


Jl47l r lJ 


U.OOO^J: 


U.OOZo 






U.04toO 


/l49(*l) 


0.005544 


0.002821 


0.004466 


0.004296 


0.004281 


hi (kb)- 1 


0.475 


0.915 


0.645 


0.585 


0.650 


Na(ti) 

N 3 m(0) 

iVu(O) 
/l47(tl) 


0.03177 
12.0 xl0~ 6 
0.5532 


0.02666 
50.5 xlO" 6 
0.5323 


0.02969 
32.1 xl0~ 6 
0.5404 


0.03041 
22.6 xl0~ 6 
0.5483 


0.02963 
27.6 xlO" 6 
0.5434 


/l49(*l) 


0.005542 


0.002836 


0.004470 


0.004289 


0.004283 


(7149 (kb) 


94 


86 


81.5 


96 


85 



in fig. HI From it we obtain two solutions for the resonance energy shift A r : 

- 11.6meV < A r < +26.0meV (30) 

- 101.9meV < A r < -79.6meV (31) 

The left-branch negative solution is tentatively ruled out by Fujii et al. j(| on the basis of 
other data, but we retain it here for completeness. 

Damour and Dyson [^j analyzed the dependence on a of the Coulomb energy contri- 
bution to the resonance energy of 149 Sm and concluded that an energy shift due to the 
electromagnetic interaction was related to a shift in a by 

A r 



Act 

a 



1.1 MeV 



(32) 



The energy 1.1 ±0.1 MeV was estimated from known isotopic shifts in the mean square 
radii of the proton density distributions in samarium. 

With our results for the bounds on the 97.3-meV resonance shift A r in 149 Sm in Eq. 
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TABLE VI: RZ2 borehole SC36 calculated fluences <f)t\ (column 4) and ancient cross sections oj 
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(column 8). Input data for Sm and U isotope fractions (columns 2,5,6 and 7 are from Ref. 
For these calculations h = 850 kyr. The metasample input data are an average of the data for the 
five individual samples. 



Sample 


N 5 (today) 
Ng (today) 


iV 6 (ti) 
JVs(*i) 


hi (kb)" 1 


f 144(h) 


fMh) 


f 149(h) 


0-149 (kb) 


1408 


0.00545 


0.02864 


0.905 


0.0046 


0.506 


0.0036 


69 


1410 


0.00527 


0.02767 


1.015 


0.0015 


0.532 


0.0030 


78 


1413 


0.00410 


0.02153 


1.820 


0.0007 


0.534 


0.0018 


65 


1416 


0.00522 


0.02741 


1.045 


0.0003 


0.534 


0.0024 


91 


1418 


0.00574 


0.03013 


0.745 


0.0003 


0.541 


0.0042 


75 


met a 


0.00516 


0.02708 


1.085 


0.0015 


0.529 


0.0030 


71.5 



we therefore obtain 2cr bounds on the a-change of: 

Aa 



0.24 x icr 7 < 



< 0.11 x TO -7 



a 



or 



0.72 x 10~ 7 < — < 0.93 x 10" 

a 



(33) 



(34) 



Our result is quite insensitive to the assumption of fixed p. Varying the H/U ratio changes 
p, and therefore affects the neutron spectrum and the resulting effective cross section. The 
major effect is due to change in the overlap of the neutron spectrum with the 149 Sm 97.3- 
meV resonance. Combining the effect on the derived 01,49 an d the change in the overlap of 
the neutron spectrum, we find S(A r ) ~ (50meV) • Sp. A 10% change in p therefore leads to 
a change in A r of order 5 meV, well within the 2a bounds assumed in extracting our final 
result. 

In their analysis Damour and Dyson assumed that changes in the nuclear part of the 
Hamiltonian did not correlate with with changes in the Coulomb energy over time. They 
noted that changes in the m q /m p could, in principle, show themselves as changes of the 
resonance position. However, at that time there was no theory to allow estimates of such 
effects. Recent QCD-based developments have speculated that a time variation in a could 
be accompanied by a larger (up to__ a factor of 30 or more) variation of the QCD scale 
parameter A (= 213 MeV) [TJ, 
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261 ] which characterizes the masses of participating 
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FIG. 4: Expanded version of the RZ10 cross section plot of fig. [3] for temperatures 200 C and 300 
C, and the bounds on A r obtained from our ±2cr limits on the (T149. 

particles. Changes in the effective nuclear potential from, e.g. changes in m 7T /m p could, in 
principle, therefore have a significant effect on the resonance shift. A detailed theoretical 
analysis and quantitative estimates of the nuclear physics aspects of the neutron resonance 
shift remain to be carried out. In such a situation the present result from the Oklo reactor 
data while stringent - and consistent with no shift in a over a two billion year period - should 
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be regarded as a model dependent 
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